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The adsorption of large ions from solution to a charged surface is investigated theoretically. 
A generalized Poisson-Boltzmann equation, which takes into account the finite size of the ions is 
presented. We obtain analytical expressions for the electrostatic potential and ion concentrations 
at the surface, leading to a modified Grahame equation. At high surface charge densities the ionic 
concentration saturates to its maximum value. Our results are in agreement with recent experiments. 
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The interaction between charged objects (interfaces, 
colloidal particles, membranes, etc) in solution is strongly 
affected by the presence of an electrolyte (salt) and is 
of great importance in biological systems and indus- 
trial applications jjJD. The main effect is screening of 
the Coulomb interaction characterized by the so-called 
Debye-Hiickel screening length || , which depends on the 
ionic strength of the solution. The Dcryaguin-Landau- 
Verwey-Overbeek theory, based on the competition be- 
tween screened Coulomb and attractive van der Waals 
interactions, has been very successful in explaining the 
stabilization of charged colloidal particles . 

One of the most widely used analytical method to 
describe electrolyte solutions is the Poisson-Boltzmann 
(PB) approach ||. For low electrostatic potentials (less 
than 25 mV), the PB equation can be linearized and 
yields the Debye-Hiickel theory ||. The PB is a con- 
tinuum mean-field like approach assuming point- like ions 
in thermodynamic equilibrium and neglecting statistical 
correlations. This theory has been successful in predict- 
ing ionic profiles close to planar and curved surfaces and 
the resulting forces. However, it is known to strongly 
overestimate ionic concentrations close to charged sur- 
faces. In particular, this shortcoming of the PB theory is 
pronounced for highly charged surfaces and multivalent 
ions. 

Since the PB equation does not take into account the 
finite size of the adsorbing ions, the ionic concentration 
close to the surface can easily exceed the maximal allowed 
coverage by orders of magnitude. Several attempts have 
been proposed to include the steric repulsion in order to 
improve upon the PB approach ||f7j. One of the first 
attempts to incorporate steric effects is the Stern layer 
modification of the PB approach. Steric effects are 
introduced by excluding the ions from the first molecular 
layer close to the surface. However, it seems difficult to 
improve on this method in a systematic way. More recent 
modifications J|f7]|] [ll]] rely either on Monte Carlo com- 



puter simulations or on numerical solutions of integral 
equations (the "hypernetted chain" equation Q). These 
approaches involve elaborate numerical calculations and 
lack the simplicity of the original PB approach. 

In this Letter, we propose a simple way to include 
steric effects in the original PB approach. This modi- 
fied PB equation clearly shows how ionic saturation takes 
place close to a charged surface. The equation is derived 
for l:z asymmetric and z:z symmetric electrolytes. At 
low ionic concentration, the original PB equation is re- 
covered. Simple analytical relations between the surface 
charge density and the counterion concentration at the 
surface are obtained, in agreement with recent experi- 
ments |||. 

Consider an asymmetric electrolyte consisting of nega- 
tive multivalent ions of charge —ze, and positive monova- 
lent ions of charge e, where e is the electron charge. The 
bulk concentration of the negative and positive ions is Cb 
and zcf,, respectively, as implied by charge neutrality. 

Within mean-field approximation, the total free en- 
ergy, F = U — TS, can be written [|l3| in terms of the 
local electrostatic potential ip(r) and the ion concentra- 
tions c (r). The electrostatic energy contribution U is: 



U 



dr 



-/i+c 1 



Stt 



\Vtp\ 2 + ec + ^p — zee ip 



(1) 



The first term is the self energy of the electric field, where 
£ is the dielectric constant of the solution. The next 
two terms are the electrostatic energies of the ions, and 
the last two terms couple the system to a bulk reservoir, 
where fj,± are the chemical potentials of the ions. 
The entropic contribution —TS is 
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where k B T is the thermal energy. For simplicity, we as- 
sume that both types of ions have the same size a. The 
first two terms are the entropies of the positive and neg- 
ative ions, whereas the last term is the entropy of the 
solvent molecules. Indeed, this last term is responsible 
for the novel steric corrections to the PB equation. In 
a more rigorous way, these corrections are obtained by 
considering a lattice-gas version of the Coulomb gas in 
which each lattice site is occupied at most by one ion 
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The variation of the free energy F = U — TS with re- 
spect to ip and yields our modified PB equation for 
the l:z electrolyte: 
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where cpQ = (z + l)a 3 Cb is the total bulk volume fraction 
of the positive and negative ions. 

For a symmetric z:z electrolyte, one gets 
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where (po — 2a c&. In the limit of small ionic concentra- 
tions, <p — > 0, Eqs. ||,[| reduce to the standard PB equa- 
tions. Moreover, for any ionic concentration and at low 
electrostatic potentials, \f3eip\ <C 1, both equations re- 
duce to the linearized PB equation (Debye-Hiickel limit) 
V 2 ip = n 2 ip where kT 1 is the Debye-Hiickel screening 
length. For the asymmetric case k 2 — 4irli,z(z + l)ch, 
where 4 = e 2 /efcsT is the Bjerrum length equal to 7 rA 
for aqueous solutions at room temperature. 

Our approach deviates significantly from the original 
PB equation for large electrostatic potentials \fleip\ ^> 1. 
In particular, the ionic concentration is unbound in the 
standard PB approach, whereas here it is always bound 
by 1/a 3 ("close packing") as can be seen from Eqs. |]]|. 
This effect is important close to strongly charged surfaces 
immersed in an electrolyte solution. 

Note that for high positive potentials, (3eip ^> 1, the 
contribution of the positive ions is negligible and the neg- 
ative ion concentration follows a distribution reminiscent 
of the Fermi-Dirac distribution [fill, 
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where the excluded volume interaction plays the role of 
the Pauli principle. 

To demonstrate the usefulness of our method, we study 
the case of a single planar surface with charge density 
a > in contact with an electrolyte solution. Ionic con- 
centration profiles are obtained from the numerical so- 
lution of Eq. U as a function of x, the distance to the 



positively charged surface. Since the positive ion con- 
centration is small near the surface, we show in Fig. la 
only the negative ion profiles, as well as the correspond- 
ing original PB profile. The main effect is the saturation 
of the ionic concentration in the vicinity of the charged 
surface. This should be contrasted with the original PB 
scheme which leads to extremely high and unphysical val- 
ues of cj = c _ (0), especially for multivalent ions. In the 
saturated region, the ionic concentration tends to 1/a 3 , 
leading to more pronounced deviations from PB for large 
ions. 

In the saturated layer the right-hand side of Eq. || 
becomes a constant, and the electrostatic potential is 
quadratic 
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where ip s is the surface potential and the boundary con- 
dition ip'\ s = —Ana/£ is satisfied. As can be seen in 
Fig. lb, the parabolic profile of ip(x) is a good approxi- 
mation close to the surface. The width of the saturated 
layer I* is not strictly equal to a. It can be easily esti- 
mated from Eq. |^ to be I* ~ a 3 cr/ze in agreement with 
Fig. la. 

The surface potential ip s can be calculated in a closed 
form from the first integral of Eq. ^, assuming that the 
concentration of the positive ions is negligible at the sur- 
face 
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Similarly, the concentration of negative ions at the sur- 
face can be calculated leading to a modified Grahame 
Equation jl| 
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This contribution is depicted in Fig. 2a, where cj is plot- 
ted as a function of the surface charge density, a/e, for 
two different ion sizes, a. The PB case is shown as well 
for comparison. At low surface charge ( < 1, the ion 
concentration reduces to the PB results 
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but for high surface charge ( > 1, the deviations from 
the PB case are substantial. Furthermore, as can be seen 
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from the above equation, the ionic concentration near 
the surface depends only weakly on the bulk electrolyte 
concentration, q,. 

In Fig. 2b the ratio between the ion charge density 
of the first layer <j\ ~ zec~a and a is plotted as func- 
tion of the specific surface area per unit charge, as is 
often measured in experiments. For the PB approach, 
this ratio diverges at high surface charge densities (small 
specific area) because c~ ~ a 2 . However, in our case, the 
steric effect changes the situation altogether since it pre- 
vents the ions from approaching and over-compensating 
the surface charges. 

The theoretical results presented here are relevant to 
recent experiments [ fp2]| where large multivalent ions are 
adsorbed onto a charged Langmuir monolayer. Large 
polyanions such as H3PW12O40 (phosphotungstic acid) 
dissolved in an aqueous subphase are attracted to a 
cationic Langmuir monolayer such as a fatty amine sur- 
factant (C20H41-NH2), spread at the water/air interface. 
The adsorbed ion density (per unit area) in the solution, 
<j\ , is measured by X-ray reflectivity. It is then related to 
the surface charge density a, which is controlled by Lang- 
muir trough. The experiments show very clearly the pres- 
ence of the steric effects for these large ions (of estimated 
size of 10 rA). As the surface charge density increases, 
<7i I a decreases in accord with our findings (Fig. 2b) and 
in contrast to the original PB approach. 

In conclusion, we have derived a modified PB equa- 
tion including steric effects. As a result the ionic concen- 
tration cannot exceed a saturation value of 1/a 3 . This 
effect is in particular important for large ions adsorb- 
ing on charged surfaces. We have obtained analytical 
expressions for the potential and ion concentrations at 
the surface. Our results differ qualitatively and quanti- 
tatively from the standard PB equation and agree with 
recent experiments on large polyions in solution. 

It would be interesting to further explore the connec- 
tion between our analytical approach and the Stern layer 
approach, as well as look at the implications on the forces 
between two charged planar surfaces. 
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FIGURE CAPTION 

Fig. 1: (a) Concentration profiles of negative multivalent ions c~(x) near a positively charged surface as obtained 
from the numerical solution of Eq. || for two different ion sizes a = 7.5 rA and a = 10 rA. Note that the saturated 
layer width is I* ~ 2 rA and 5 rA, respectively. The solid line represents the concentration profile of the standard 
PB equation, (b) Calculated electrostatic potential profiles near the surface plotted together with the parabolic 
approximation (Eqs. ||, The dotted, dashed and solid lines are as in (a). The bulk concentration is Cf, = 0.1M for 
a l:z electrolyte with 2=4. The surface charge density a is taken as one electron charge per 50 rA 2 . The aqueous 
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solution with s = 80 is at room temperature. 

Fig. 2: (a) Surface concentration of ions as a function of the surface charge density for different ions size, a. The 
PB concentration is also plotted for comparison, (b) Ratio of the first layer charge density and the surface charge 
density, <j\/o — zecja/a, as a function of the specific surface area per unit charge, e/cr, for different ion sizes. The 
PB result is plotted with the first layer width taken as 5 rA. The l:z electrolyte bulk concentration is q, = ImM and 
the valency z = 4. 
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